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We report on numerical studies of the energy spectrum and the localization properties in the 
double-layer quantum Hall system at v = 1. The Coulomb interaction is treated by the Hartree- 
Fock approximation, and the localization properties in the presence of disorder are studied by 
evaluating participation ratios for the Hartree-Fock eigenfunctions. We show that the extended 
states seem to exist only near each center of the two subbands split by the exchange-enhanced 
energy gap. It is also shown that the self-consistent orbitals whose energies are close to the Fermi 
energy appear to become extended together with the reduction in the energy gap as the layer 
separation increases. The collapse of the energy gap expected from our results is consistent with 
the incompressible-compressible transition observed in recent experiments, and the change of the 
localization properties near the Fermi energy can explain the disappearance of the quantum Hall 
effect for large layer separations very well. 



I. INTRODUCTION 

When the integer quantum Hall effect (QHE) is stud- 
ied theoretically, it is usually assumed without justifi- 
cation that the Coulomb interaction between electrons 
can be safely ignored. In a strong magnetic field, the 
eigenfunctions of the single-particle Hamiltonian for two- 
dimensional non-interacting electron systems are local- 
ized by a disorder potential at almost all energies except 
for a discrete set of critical energies {e c n} near the cen- 
ter of each disorder-broadened Landau level. Theoretical 
studies suggest that at T = the Hall conductivity jumps 
by e 2 /h each time the Fermi energy sp crosses one of the 
critical energies, and that the longitudinal conductivity 
is zero if Sf £ c n- These suggestions are supported by 
many experimental studies jl| . 

However, such localization properties are not always 
guaranteed when two different Landau levels are nearly 
degenerate. In fact, the numerical studies of the double- 
layer quantum Hall (QH) system in the absence of 
electron-electron interactions could not obtain the 
reasonable localization properties in case of nearly degen- 
erate Landau levels (i.e., in the weak interlayer-tunncling 
case). This is a typical case where the Coulomb inter- 
action should be considered even for understanding the 
integer QHE qualitatively. Thus we consider this system 
in the presence of interactions in this paper. We show 
that the exchange-enhanced energy gap Q appears in 
this interacting system and that the localization prop- 
erties consistent with the observation of the QHE can 
be obtained for small layer separations even in the weak 
interlayer-tunneling case. 

In double-layer QH systems ||, the interlayer- 
tunncling of electrons brings about the mixing of the Lan- 
dau levels in the two layers, and the Landau levels split 
into symmetric and antisymmetric combinations about 



the center of the double-layer structure. The energy gap, 
AsaSj between them is proportional to the tunneling am- 
plitude, and it needs to be small for the nearby degener- 
acy of the Landau levels. 

Such samples are realized experimentally and the 
transport properties have been investigated HJt]]. In 
these experiments using high-mobility samples, it has 
been reported that the Coulomb interaction plays an im- 
portant role on the ground-state properties and low-lying 
excitations. For v = 1, the phase diagram against the 
layer separation, d, which controls the strength of the 
interlayer interactions, and Asas is obtained experimen- 
tally 0. The phase diagram shows that the QH state 
disappears for d > d c and that the critical separation d c 
increases as Asas increases. 

In a more recent experiment using high mobility sam- 
ples with weak interlayer-tunneling, the zero in the lon- 
gitudinal resistivity is replaced for d > d c by a broad 
minimum similar to that observed in the single-layer QH 
system at v = 1/2 ||. This suggests a transition from 
an incompressible QH state with strong interlayer cor- 
relations to compressible state consisting of two (weakly 
correlated) layers, where the metallic states of composite 
fermions are formed. 

Theoretically the pseudospin formalism is often in- 
troduced to describe the layer degrees of freedom in 
double-layer systems. This is done by assigning the up- 
per/lower layers to the pseudospin T / 4- At f = 1/m 
(to: odd integer), the pseudospin ferromagnetism re- 
sults from the interlayer-tunncling and exchange interac- 
tions between electrons . The phase boundary between 
the QHE and the non-QHE phase was determined the- 
oretically by assuming that the QHE phase is destroyed 
together with the collapse of the pseudospin ferromag- 
netism P)^o[. The pseudospin- ferromagnetic ground 
state is shown to evolve continuously from tunncling- 
dominated to correlation-dominated as Asas decreases 
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The Hartree-Fock calculations have also been done to 
study the QH systems, and this approximation is ex- 
pected to describe the electronic properties well espe- 
cially for the integer filling factors. In fact, this approx- 
imation has been used for the study of the double-layer 
v = 1 QH system in the absence of random disorder po- 
tential 

In this paper, we investigate the energy spectrum and 
the localization properties in disordered double-layer QH 
systems at v = 1. The Coulomb interaction is treated 
by the Hartree-Fock approximation, and the localization 
properties are studied by evaluating participation ratios 
for the Hartree-Fock eigenfunctions. This method has 
been used in the studies on the interaction effects in the 
single-layer QH system [|l4|Jl5H . We show that the lo- 
calization properties change qualitatively because of the 
interaction effects JlfJ . 

Our paper is organized as follows. In Section II, we 
explain the model and calculation methods which we 
use in this study. In Section III, we first discuss the 
results of the numerical calculations in the absence of 
Coulomb interactions. Through this discussion about 
the previous results m and ours, it is shown that 
the Coulomb interaction should be considered to under- 
stand the v = 1 QHE in double-layer system with weak 
interlayer-tunneling. After this discussion, we show our 
numerical results in the presence of Coulomb interac- 
tions. These are our main results in this paper, and it 
is shown that the localization properties change together 
with the reduction in the exchange-enhanced energy gap 
as the layer separation increases. Finally in Section IV, 
we briefly summarize our findings. 



In the strong- magnetic-field limit, it is enough to con- 
sider only the lowest Landau level because one can ne- 
glect the Landau level mixing by interactions or disor- 
ders. The real spin degrees of freedom is also ignored 
by assuming the spin polarization due to the Zeeman 
energy. Our attention is restricted to this strong-field 
limit. We apply the periodic boundary condition to the 
single-particle wavefunctions inside the two parallel rect- 
angles of dimensions L x , L y , and use the Landau gauge 
A(r) = (0, Bx, 0). One can then use the following set of 
basis functions for the lowest Landau level: 
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where I = yJch/eB is the magnetic length, and Xj — 
2nt 2 i [L y is the center coordinate of the j-th Landau or- 
bit WJ\. The orbital de grees of freedom in each layer 
are described by this set of the Landau orbits, and the 
two layers can be distinguished by the pseudospin index 
<7. Thus the set {\ja)} (j = 1, 2, • • • , N L , a = f, j) can 
be used as a basis set for representing the Hartree-Fock 
Hamiltonian, where Nl — L x L v /2tt£ 2 is the Landau level 
degeneracy in each layer. 

In terms of the set of basis functions, {|jcr)}, the ma- 
trix element of the Hartree-Fock Hamiltonian is given by 



(ja\HuF\f<r') = 7T^"%' + <W (jo\Vx mp \j'o-) 



(1) 



+0v|vw|jV), 



(2) 



II. MODEL AND METHOD 

We consider a double-layer system of spin-polarized 
electrons in a strong magnetic field perpendicular to 
the layers. In double-layer systems, there exists the 
interlayer-tunneling of electrons. The single-particle 
wavefunctions then split into symmetric and antisymmet- 
ric ones about the center of the double-layer structure, 
and the energy gap between them, AgAS, enters as an 
energy scale. The thickness of the wavefunction in each 
layer is neglected for simplicity. 

The two-dimensional coordinates in the two parallel 
planes are denoted by r = (x,y), and the layer degrees 
of freedom are described by the pseudospin a = t> I- 
The Coulomb interaction between electrons is then de- 
pendent on pseudospin a for a finite layer separation d. 
Its Fourier transform V aa i (q) is 2we 2 /eq if a = a' (i.e., for 
the intra-layer interaction) and (2ire 2 /eq) e~ qd if a ^ a' 
(i.e., for the inter-layer one), where e is the dielectric 
constant of the host material. The Coulomb interaction 
is treated self-consistently within the Hartree-Fock (HF) 
approximation. 



where Sjj> is a usual Kronecker delta, the first term in 
the right hand side of eq.(||) is due to the interlayer- 
tunneling. The second term and the third one in eq.(Q) 
result from the impurity scattering and the Coulomb 
interaction between electrons, respectively. The ampli- 
tudes of interlayer impurity scatterings are neglected be- 
cause of their small values and for simplicity. 

The Hartree-Fock single-particle equation is given by 

#HF |<y9a) = £ a \<p a ) , (3) 

where e a and \<p a ) are an eigenvalue and correspond- 
ing eigenstate of this equation, respectively. Because the 
third term, (j&\VHF\j'a'), in eq.(g) is dependent on the 
set, {|y a )}, as seen in eq.([|), this single-particle equation 
must be solved self-consistently. This is done by diago- 
nalizing the 2Nl X 2Nl matrix, (ja\H-HF\j'a'), in eq.(g) 
numerically and solving eq.(|J) iteratively until the self- 
consistency is achieved. Among the obtained eigenstates, 
{|y a )}j m eq-©! the N lowest-energy ones are occupied 
in A^-electron systems, and N = Nl in case of v = 1. 

The matrix element of the impurity scattering, which 
is the second term in eq.(||), is given by 
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where v a (q) is the Fourier transform of the impurity po- 
tential, v a (r), in the layer a, S(j, j') is 1 if j = j' (mod 
Nl) and otherwise. The sum over the wave vector q 
is over q x = (2ir/L x ) n x , q y = (2ir/L y )n y (n x , n y : inte- 
gers), because the periodic boundary condition is used. 

Our model disorder consists of randomly located 5- 
function scatterers with a random strength uniformly dis- 
tributed between —Vq and Vq. The disorder potential is 
then given by 



Wtr (r) = £^(r-Rn, 



(5) 



where Vf and Rf are the strength and position of the 
j-th impurity in the layer a, respectively. There exist 
iVi m p impurities in each layer, and we assume that the 
disorder potentials in the two layers are uncorrelated, i.e., 
there is no correlation about {Rf } and {V?} between the 
two layers ||. For this model disorder, the energy scale 
which characterizes the Landau subband width is given 
by T = (ViN imp /e 2 L x L y )^ 2 @. We choose to work 
with N imp /N L = 2tt£ 2 N imp /L x L y = 5, and keep V /£ 2 
constant in order to use T as the unit of energy. 

The matrix element of the Coulomb interaction, 
(j^\VnF\j'cr'), which is the third term in eq.(|2|), is given 
by 
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where the sum over q in cq.(|^) is over 'the Brillouin 
zone', q x = (2n/L x )n x , q y = (2ir/L y )n y (n x , n y = 
1,2,---,N L ), and the quantities Z7£ CT '(q) and Ug a ' {q) 
correspond to the Hartree and Fock potential, respec- 
tively [|l9) . The quantity A aa > (q) is proportional to the 
expectation value of the density operator p aa >(q). The 
two-body matrix element, (jihlV^'ljaji), in eq.(jl^) is 
determined by the Fourier transform, V aa '(q), of the 
Coulomb interaction, and a field operator, ^ a (r), for 
pseudospin a is considered within the subspace of the 
lowest Landau level. In eq.(^|) and eq.^), 6{x) and ep 
are the Heaviside step function and the Fermi energy, re- 
spectively, and only the N lowest-energy eigenstates in 
eq.(||) contribute to the sum over a. About the quantities 
Up a (q), (q), and A CTCT '(q), it is enough to consider 
them only within 'the Brillouin zone' because of their 
periodicity. 

When the HF single-particle equation is solved, 
the quantity A aa > (q) can be used to check the self- 
consistency of the calculated results. We have judged 
the convergence of the calculated ones by the following 
condition 



AN? <H ^ 

L cr,a' qGB.Z. 



(12) 



where k represents each iteration step. We note that our 
calculations have been done under the constraint that 
the average numbers of electrons are the same in the two 
layers, i.e., A CT(T (0) = v/2 = 1/2 for a = T, J. |1- 

In order to investigate the localization properties, we 
evaluate the participation ratios for the self-consistent 
Hartree- Fock eigenstates [jl5|. The participation ratio is 
given by 



P a = 



L X Ly 



dx / dy \ip a (r) 



(13) 



for a normalized eigenstate y a (r) = (r\ip a ), and P a ~ 
^/L x L y , where £, a is the localization length of an eigen- 
state, y> a (r). As the eigenstate becomes more extended, 
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the participation ratio becomes larger. That is, the par- 
ticipation ratio shows how extended the eigenstate is. 

For our numerical calculations, Asas and e 2 /e£ can 
be used as the energy scales for the interlayer-tunncling 
and the Coulomb interaction, respectively. As the unit 
of energy, the strength, T, of the impurity potential is 
used, and we consider the following cases, AsAs/r = 0.1, 
(e 2 /el)/T = 20, d/£ = 1.2, 1.5, 1.8. For these values of the 
parameters, the amplitude of the interlayer-tunneling is 
very small and the Coulomb interaction is much stronger 
than the disorder potential (e 2 /e£ ~ 100K for typical 
GaAs samples) . Thus we can compare the calculated re- 
sults with the experimental ones for high mobility sam- 
ples with weak interlayer-tunncling. 

In the process of our self-consistent calculations, we 
first solve eq.(^) in the absence of Coulomb interactions, 
i.e., for (e 2 /e£)/r = 0, where no self-consistency is re- 
quired. By increasing the ratio, (e 2 /e£)/T, gradually to 
the required value and using the latest results as input 
data in a new calculation step JEj] , we have obtained the 
self-consistent solution of eq.(^)~^ 

We also show the results in the absence of Coulomb 
interactions, because the comparison between them and 
those in the presence of interactions clarifies the impor- 
tance of Coulomb interactions for the localization prop- 
erties in the double-layer QH system at v = 1. The 
calculations in the absence of interactions are done for 
AsAs/r = 1 and 0.1, i.e., in both strong and weak tun- 
neling cases. 



III. NUMERICAL RESULT AND DISCUSSION 
A. In the absence of Coulomb interactions 

We first discuss the double-layer QH system in the ab- 
sence of Coulomb interactions. The localization prop- 
erties in this non-interacting system were studied previ- 
ously ^,0]. Although our numerical results in the non- 
interacting case are almost similar to previous ones j|, 
we show them in Figure 1 to make our discussion easy to 
understand. The importance of Coulomb interactions in 
double-layer QH systems can be understood by compar- 
ing these results with ones in the presence of interactions, 
which will be given in the next subsection. 

We calculated the density of states (DOS) and partic- 
ipation ratio for several impurity configurations in the 
non-interacting double-layer QH system with the Lan- 
dau level degeneracy, Nl = 256. Because the square 
(L x = L y = L) systems are now considered, the dimen- 
sion L is given by L — — 40£ {£: the magnetic 
length) . Figure 1 shows the result for one of these impu- 
rity configurations in the cases of AsAs/r = 1 and 0.1. 
In the non-interacting case, the parameter, AgAs/r, is 
needed to characterize the disordered double-layer QH 
system. The layer separation, which controls interlayer 



interactions, does not need to be specified in spite of its 
importance in the real experiments. 

In the figures of the DOS in Figure 1, the numbers 
of eigenstates within a finite width of energy, AE, are 
plotted. The width AE is several times as large as the 
average energy-level spacing. We note that the energy E 
in the figures is given in units of T all through this paper 
and that each of the vertical broken lines indicates the 
highest energy eigenvalue among those of the occupied 
eigenstates. 

For AgAs/r = 1, which corresponds to the strong 
interlayer-tunneling case, there exist two subbands which 
mainly consist of symmetric and antisymmetric combina- 
tions of isolated layer states, respectively. The localiza- 
tion properties within each subband are similar to those 
in the single-layer QH systems as pointed out previously 
That is, the eigenfunctions are extended only near 
the center of each disorder-broadened subband. 

For AgAs/r = 0.1, however, the symmetric and anti- 
symmetric subbands are not well developed in the DOS, 
and the localization properties are quite different from 
those for AgAs/r = 1. In fact, the peak in the par- 
ticipation ratios has a much broader width than that 
for AgAs/r = 1, and the participation ratios near the 
Fermi energy take much larger values than those for 

AgAs/r = i. 

In the previous study || , the finite-size scaling method 
was also used for the non-interacting system with weak 
interlayer tunneling. Then it was claimed that the ex- 
tended states exist only at the two energies, which are 
split by somewhat larger than AgAs, rather than across 
a band of finite width between the low- and high-energy 
mobility edges in the thermodynamic limit. In the en- 
ergy interval between the two extended-state energies, 
however, their numerical values of the localization length 
in the weak-tunneling case are much larger than those 
in the strong-tunneling case. Moreover it is not clear 
from the system-size dependence of their results whether 
the extended states do not exist across a band of finite 
width in the thermodynamic limit, either. Thus their 
claim seems to be controversial as far as one judges from 
their numerical results. Therefore we consider the effects 
of the Coulomb interaction ignored in their study, and 
show that the extended states seem to exist only near 
each center of the two subbands split by the interaction. 



B. In the presence of Coulomb interactions 

In the presence of Coulomb interactions, it is necessary 
to specify the layer separation, d/£, which controls inter- 
layer interactions, as well as the interlayer-tunneling am- 
plitude, AgAs/r, in order to characterize the disordered 
double-layer QH system. In this subsection, we consider 
the double-layer QH system with AgAs/r = 0.1, i.e., the 
weak interlayer-tunncling case. The localization proper- 
ties in this case remain unclear as seen in the previous 
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subsection. About the layer separation, we consider the 
three cases of d/i = 1.2, 1.5, 1.8. By using these values, 
we can see that the localization properties change quali- 
tatively as the layer separation increases. 

We calculated the DOS and participation ratio for sev- 
eral impurity configurations in the double-layer QH sys- 
tem with the Landau level degeneracy, Nl = 80. The di- 
mension L is then given by L = £^/2itNl — 22i. Figure 2 
shows the result for one of these impurity configurations 
in the cases of A SAS /r = 0.1, d/i = 1.2, 1.5, 1.8. In the 
figures of the DOS, the numbers of eigenstates within a 
finite width of energy, AE, are plotted, and the energy 
E in the figures is given in units of T. Each of the ver- 
tical broken lines in the figures of the participation ratio 
indicates the highest energy eigenvalue among those of 
the occupied eigenstates. 

Let us first consider the case of d/£ = 1.2 in Figure 2. 
We can see that the DOS has a large energy gap near 
the Fermi energy. The participation ratios take small 
values (nearly zero) at the edges of each subband, and 
take much larger values around the center of each sub- 
band. These localization properties are almost similar to 
those for AsAs/r = 1 in Figure 1, and are quite different 
from those for AsAs/r = 0.1 in Figure 1. This difference 
results from the electron correlation effects and the en- 
ergy gap in Figure 2 is due to the exchange interactions 
between electrons. In fact, for small layer separations, 
d/i , the ground state is a pseudospin-ferromagnetic one 
jfl| , and the symmetric and antisymmetric combinations 
of isolated layer states are separated from each other by 
an exchange-enhanced energy gap. Then the localiza- 
tion properties within each subband are expected to be 
similar to those in the single-layer QH systems and be 
consistent with the observation of the integer QHE. 

For larger layer separations [dji — 1.5, 1.8), however, 
the localization properties become complicated. The par- 
ticipation ratios for the electronic eigenstates whose en- 
ergies are close to the Fermi energy take larger values 
than those in case of d/i = 1.2, i.e., they are more ex- 
tended than the ones for d/i = 1.2. Moreover, the energy 
gap seen in the DOS seems to decrease as the layer sep- 
aration increases. This decrease in the energy gap corre- 
sponds to the theoretical results from the viewpoint of the 
pseudospin ferromagnetism PJlOfl that the pseudospin- 
ferromagnetic order are broken gradually with the in- 
crease of layer separation. If the energy gap continues 
to decrease and eventually collapses with the increase of 
layer separation, it is consistent with the incompressible- 
compressible transition reported in recent experiments 

Although the size-scaling calculations are needed to 
make the localization properties clear, we claim the exis- 
tence of the extended states only near the center of each 
subband in case of d/i = 1.2 from the following facts. 
One of them is that the dependence of such localization 
properties on the layer separation is consistent with the 
transition observed experimentally between the QHE and 
the non-QHE phase @,^|. The other is that the localiza- 



tion properties in case of d/i — 1.2 are almost similar to 
those in case of Asas /r = 1 in Figure 1, although the ori- 
gins of the energy gap are different. We expect that our 
claim will be confirmed by the size-scaling calculations. 

We also performed numerical calculations in the cases 
of A S As/r = 1, d/i = 1.5,1.8,2.0,2.3. Although 
these results are not shown graphically in this paper, 
we obtained the results which are similar to those for 
AsAs/r = 0.1. The energy gap then survives for larger 
layer separations than that in case of AsAs/r = 0.1. 
This corresponds to the experimental and theoretical 
results [^,0 that the critical separation d c increases as 
Asas increases. 

Thus the exchange-enhanced energy gap under the 
pseudospin-ferromagnetic order appears between the two 
subbands, and the existence of the extended states only 
near the center of each subband seems to be realized 
for small layer separations by the electron correlation 
effects. Unfortunately our Hartree-Fock calculations in 
this paper are limited to those for small system sizes and 
relatively small layer separations, because it is difficult 
to obtain the numerical convergence for larger system 
sizes or larger layer separations. Therefore, for the con- 
firmation of our claim about the localization properties in 
weak interlayer-tunneling case, further numerical studies 
are needed for larger systems. Especially the size-scaling 
calculations should be done to discuss the change of the 
localization properties quantitatively. 



IV. SUMMARY 



We investigated the disordered double-layer QH sys- 
tem at v — 1 numerically, and then the Coulomb inter- 
action was taken into consideration within the Hartree- 
Fock approximation. We examined the density of states 
and the participation ratios for the self-consistent eigen- 
states of the Hartree-Fock Hamiltonian. By considering 
the interaction effects in this disordered system, it was 
found that the extended states seem to exist only near 
each center of the two subbands split by the exchange- 
enhanced energy gap for small layer separations. It was 
shown that the eigenstates whose energies are close to 
the Fermi energy appears to become extended together 
with the decrease in the energy gap as the layer separa- 
tion increases. The collapse of the energy gap, which can 
be expected from our results, is also consistent with the 
incompressible-compressible transition reported in recent 
experiments. Especially the change of the localization 
properties near the Fermi energy can explain the disap- 
pearance of the QHE phase for large layer separations 
very well. 
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FIG. 1. The density of states (DOS) and participation ratio in the non-interacting double-layer QH system at v = 1 are 
shown in both strong (AsAs/r = 1) and weak (Asas/F = 0.1) interlayer-tunneling cases. In the figures of the DOS, the 
numbers of eigenstates within a finite width of energy, AE, are plotted. Each of the vertical broken lines in the figures of 
the participation ratio indicates the highest energy eigenvalue among those of the occupied eigenstates. The Landau level 
degeneracy, Nl, is 256, and the energy E in the figures is given in units of V. 
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FIG. 2. The density of states (DOS) and participation ratio for the self-consistent eigenstates in the interacting double-layer 
QH system at v — 1 are shown in the weak interlayer-tunneling case (AsAs/r = 0.1). As the layer separation increases, the 
energy gap decreases and the electronic eigenstates whose energies are close to the Fermi energy seem to become extended. 
The values used for the layer separation, d/l, are 1.2, 1.5 and 1.8, and the Landau level degeneracy, Nl, is 80. 
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